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Abstract: The complex interactions in desert ecosystems between functional types and environmental 
conditions could be reflected by plant water use patterns. However, the mechanisms underlying the water 
use patterns as well as the water sources of Tamarix laxa in the mega-dunes of the Badain Jaran Desert, 
China, remain unclear. This study investigated the water sources and water use patterns of T. axa using 
the stable oxygen isotope method. The 6!8O values of xylem water, soil water in different layers (0-200 
cm), rainwater, snow water, lake water, atmospheric water vapor, condensate water, and groundwater were 
measured. The sources of water used by T. axa were determined using the IsoSource model. The results 
indicate that T. /axa mainly relies on soil water. At the beginning of the growing season (in May), the 
species is primarily dependent on water from the middle soil layer (60-120 cm) and deep soil layer 
(120-200 cm). However, it mainly absorbs water from the shallow soil layer (0—60 cm) as the rainy season 
commences. In September, water use of T. /axa reverts to the deep soil layer (120—200 cm). The water use 
patterns of T. /axa are closely linked with heavy precipitation events and soil water content. These findings 
reveal the drought resistance mechanisms of T. /axa and are of significance for screening species for 
ecological restoration. 
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1 Introduction 


Tamarix laxa, an important shrub with the height of 1.5-3.0 m in arid and semi-arid regions, is 
widely distributed on desertified and salinized land (Wu et al., 1990). Globally, there are more 
than 90 species in the genus Tamarix, which are mainly distributed in western Europe, Asia, and 
southwestern Africa. In China, T. laxa is widely distributed in the northern part, including the 
Hexi Corridor, Tengger Desert and Badain Jaran Desert. T. laxa has a well-developed root system, 
a strong germination ability and also a high ecological adaptability, and possesses adaptations for 
wind defense, sand fixing, water and soil retention, and feed value (Yin, 1995; Horton and Clark, 
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2001; Wu et al., 2009; Polacik and Maricle, 2013). Thus, this species has a great potential for 
desertification control and ecological restoration in arid and semi-arid regions. 

T. laxa is a shrub with excellent wind and sand fixation functions, growing in saline sandy 
lands, lakeshores, canal sides, river terraces, and river floodplains. It is well known that in desert 
environments, water is the main factor limiting the distribution, growth, and development of 
vegetation. Surprisingly, we discovered that in the southeast of the Badain Jaran Desert, T. laxa 
plants with the height approximately up to 3.0-4.0 m and a ground diameter of up to 6.3 cm are 
distributed in the upper and middle sections of 14 mega-dunes near the Dashazao Lake, Hulai 
Jaran and Zhalate Lake. These plants are generally located on the leeward slopes. T. laxa is 
assumed to be a typical groundwater-dependent plant species, with a well-developed root system 
relying on groundwater to survive (Gries et al., 2003; Yu et al., 2017). However, the mega-dunes 
of the Badain Jaran Desert where T. laxa distributed have a relative height of 101—244 m, making 
access to groundwater difficult (Yang et al., 2014). Thus, it is interesting to explore the sources of 
water used by T. laxa in this region, and to understand water use patterns of T. laxa. 

Previous studies on Tamarix have mainly focused on their physiological and ecological 
characteristics, community structures and leaf functions (Peng et al., 2010; Ma et al., 2013; Zhao 
et al., 2017). The water sources of Tamarix have been explored in the Heihe River Basin, Junggar 
Basin, Qaidam Basin and Dunhuang oasis by many researches (Zhao et al., 2008; Zhou et al., 
2011; Yin et al., 2012; Zhou et al., 2013; Xing et al., 2014; Cui et al., 2017). To date, no water 
isotope study has been conducted to explore the water sources and water use patterns of T. laxa in 
desert environments, especially for mega-dunes. As a consequence of global climate change, 
significant groundwater depletion has occurred in the transition zone between deserts and oases, 
resulting in vegetation losses (Xu, 2008; Qin and Stocker, 2014). Thus, investigation into water 
sources and water use patterns of desert plants is essential to ecosystem protection and 
restoration. 

Generally, isotopic signatures for both 180/160 (8!8O) and deuterium/hydrogen ratio (5D) vary 
in different types of water resources (Dansgaard, 1964). Xylem water can reflect the isotopic 
composition of the water used by plant species (Dawson, 1993; Dawson and Ehleringer, 1998). 
Hence, stable isotope analysis is considered as a powerful method to investigate plant-water 
relationships (Flanagan et al., 1992; Dawson, 1996). This method has been widely applied in 
many fields, such as forestry, ecology and hydrology (Gazis and Feng, 2004). Researches have 
shown that in arid and semi-arid regions, water sources of plants differ among plant species, water 
environments and growing seasons (Cao et al., 2002; Dai et al., 2015; Wu et al., 2016). Moreover, 
the atmospheric water vapor may also be an important water source for desert plant species (Xiao 
et al., 2014). 

The objective of this study was to investigate the water sources and water use patterns of T. 
laxa on the leeward slopes of mega-dunes in the Badain Jaran Desert. It should be noted that we 
used oxygen isotopes instead of hydrogen isotopes to determine the water sources of T. laxa 
species, because some previous studies have suggested that the fractionation of hydrogen isotopes 
occurs as water enters the root tissue (Sternberg et al., 1991; Lin and Sternberg, 1992; Ellsworth 
and Williams, 2007). The results of this study can reveal the drought resistance mechanisms of T. 
laxa in desert ecosystems and have a great significance for ecological restoration in the Badain 
Jaran Desert as well as some other similar regions in the world. 


2 Materials and methods 


2.1 Study area 


This study was conducted in the southeast of the Badain Jaran Desert (39°20'-42°00'N, 
99°48’-104°14’E; Fig. 1), which is located in the northwestern Alagxa Plateau of western Inner 
Mongolia, China. The area of the Badain Jaran Desert is 5.2x104 km? (Zhu et al., 2010). A number 
of mega-dunes are distributed in the southeastern part of the desert, with relative heights ranging 
from 200 to 300 m (Dong et al., 2004). Some springs and groundwater exist near the desert lakes, 
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with a total dissolved solid content of less than 1 g/L. The region is dominated by an arid 
continental climate, with mean winter and summer temperatures of —9.1°C and 25.3°C, 
respectively. Mean annual precipitation is approximately 76.9 mm, increasing from 35.2—-42.9 
mm in the north to 90.1—115.4 mm in the south (Ma et al., 2011). Potential evaporation exceeds 
2000.0 mm, which is 20 times of the mean annual precipitation (Dong et al., 2016). Annual mean 
wind speed increases northward from 2.8 to 4.6 m/s, with the strongest winds occurring in April 
and May. Soil bulk density is 1.61 g/cm*. Plant species richness is low, with the dominate species 
of shrubs (Artemisia arenaria, Calligonum mongolicum, Nitraria tangutorum and Sarcozygium 
xanthoxylon) and herbs (Psammochloa villosa, Corispermum patelliforme and Agriophyllum 
squarrosum). 
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Fig. 1 Overview of the Badain Jaran Desert and locations of the four sampling plots in the Badain Jaran Desert. 
Pı, Dashazao Lake; P2, Hulai Jaran; P3, Zhalate Lake; P4, Chaxingtaoletu. 


2.2 Sampling and data collection 


This study was conducted in 2015. Four mega-dunes (relative height of about 200 m) with T. laxa 
growing on the leeward slopes were chosen as sampling plots (Fig. 1; Table 1). T. laxa has a 
well-developed root system (Fig. 2). We found that the vertical root length of T. laxa can reach up 
to 130-180 cm, and the root biomass is concentrated at the 30-120 cm soil depths. At greater 
depths, the roots are rotting and have no ability to absorb water. 


Table 1 Geographical information of the four sampling plots 


Plot name Plot number Geographical coordinate Altitude (m) Slope (°) 
Dashazao Lake P: 39°46'15"N 102°08'20"E 1344 20-23 
Hulai Jaran P2 39°36'04"N 102°02'32"E 1409 18-15 
Zhalate Lake P; 39°48'54"N 102°23'36"E 1398 25-27 
Chaxingtaoletu Py 39°46'28"N 102°23'27"E 1436 23-25 


The rainfall sample was collected using a rain gauge, which is constructed with a 10-cm 
diameter funnel connecting to a measuring bottle. Rainfall samples were collected immediately as 
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soon as the rainfall events ended to avoid evaporation. Totally, 18 event-based rainfall samples 
were collected from May to September, 2015 at Dashazao Lake (Pı) and Chaxingtaoletu (P4). 
Twenty-three atmospheric water vapor samples were collected using aluminum foil gas bags with 
volume of 1 L from May to September, 2015 around the sampling plots. The specific collection 
process was as follows: (1) closing the valve and connecting the rubber ball; (2) opening the valve 
and pinching the ball; (3) collecting the water vapor samples; (4) closing the valve; and (5) 
removing the rubber ball. Snow samples were collected on 17 January, 2015 at Alagxa Right 
Banner and sealed in plastic bottles. They were kept at room temperature until completely melted 
and then transferred into clean polyethylene bottles (50 mL). Moreover, 14 groundwater and 33 
lake water samples were respectively collected from nearby wells and lakes using polyethylene 
bottles. Six atmospheric condensation water samples were collected using a homemade 
equipment (Patent No.: ZL201620375043.2). All these water samples were sealed with parafilm 
and stored in a fridge at 4°C. 


Re) 


Fig. 2 Tamarix laxa growing in the Badain Jaran Desert (a), vertical root distribution (b) and rot roots (c) 


The growth characteristics of T. laxa in the four mega-dunes are provided in Table 2. Plant and 
soil samples were simultaneously collected during 26—29 May, 10-12 July, 20-22 August and 
21-24 September, 2015. It should be noted that no precipitation events occurred during the 
sampling periods. On each sampling date, we chose three sampling plants in each plot to collect 
suberized twigs (diameter of 6-10 mm and length of 50-60 mm). As soon as sampling, the green 
tissues, including all leaves and green bark, were removed quickly from the stems to avoid 
isotopic fractionation of the xylem water. Then, they were immediately put into screw-cap glass 
vials (volume of 8 mL), sealed with parafilm, and stored in a freezer at —20°C until water 
extraction (Dawson, 1996; Liu et al., 2014). Soil samples were collected using a hand auger close 
to the sampled plants at the same time the suberized twigs were sampled. The soil sampling depth 
was 0-200 cm and the sampling interval was 20 cm. Four replicates were collected from each soil 
layer and then they were mixed. Soil sample for each layer was divided into two portions. One 
was sealed in a glass vial and frozen with the plant sample for isotope analysis, and the other was 
stored in an aluminum box for soil water content measurement by the oven drying method. 
Cryogenic vacuum distillation method was used to extract water from plant and soil samples. Soil 
particle size fractions were measured using a laser particle sizer (Mastersizer 2000, Malvern 
Instruments, Worcestershire, UK). Soil particle size composition shows no significant variations 
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with depth (P>0.05). Soil is mainly composed of fine sand (50.63%) and medium sand (44.48%) 
according to the USDA (United States Department of Agriculture) classification (Table 3). 
Meteorological variables, such as daily precipitation, were measured by an automatic weather 


station (Vaisala, MAWS301, USA) located at Cherigele Lake (39°53'00"N, 102°15’40”E). 


Table 2 Growth characteristics of Tamarix laxa in the four sampling plots 


Ground diameter (cm) 


Plot number Height (m) Crown area (mxm) 
Maximum Minimum Average 
Pi 3:2 9.0x9.5 5.65 3.45 4.58 
P2 3.5 12.0x15.0 6.57 3.53 5.16 
P3 35) 8.0x10.0 4.85 3.20 4.20 
Py 3:3 12.5x10.8 6.91 4.15 5.82 


Table 3 Composition of soil particle size at different depths 


Silt Very fine sand Fine sand Medium sand Coarse sand 
Soil depth (cm) (0.005—0.05 mm) (0.05—0.10 mm) (0.10-0.25 mm) (0.25-0.50 mm) (0.50-1.00 mm) 
(%) 

20 0.00 3.54 60.56 34.56 1.34 
40 0.20 2.21 54.86 40.55 2.18 
60 0.37 2.21 51.73 44.84 0.85 
80 0.00 1.29 52.16 44.87 1.69 
100 0.00 2.91 48.83 45.64 2.61 
120 0.21 2.48 43.79 48.95 4.57 
140 0.00 2.07 48.84 45.17 3.91 
160 0.00 1.68 46.03 49.33 2.96 
180 0.00 2.04 53.38 41.71 2.88 
200 0.00 1.51 46.17 49.22 3.10 
Mean 0.08 2.19 50.63 44.48 2.61 


2.3 Date analysis 


The 8D and 5!80 of the extracted water from all plant and soil samples were analyzed by a water 
vapor isotope analyzer (PicarroL2120-I, USA) at the State Key Laboratory of Desertification and 
Aeolian Sand Disaster Combating, Gansu Desert Control Research Institute, China. All 
precipitation and other water samples were filtered using 0.22-um pore size filters prior to 
measurements. Isotope measurements were calibrated using the international standard, Vienna 
Standard Mean Ocean Water (V-SMOW). The final measurement accuracies for 6D and 8!8O 
were 0.1%0 and 0.5%o, respectively. The oxygen isotopic composition can be expressed as 
Equation 1: 


51 O (%0) -Zeme 1)x1000% (1) 
standard 

where Rsample and Rstandara are the oxygen isotopic compositions ('8O/'°O molar ratio) of the 

sample and standard water (Standard Mean Ocean Water, SMOW), respectively. 

The isotopic values of xylem water were compared with those of potential water sources using 
the IsoSource model to obtain a feasible range of different water sources used by T. laxa at each 
sampling time (Phillips and Gregg, 2003). We identified three potential soil water sources (i.e., 
layers of 0-60, 60-120 and 120-200 cm depth) based on the similarities of 5'°O value for each 
soil depth. Three soil layers were classified as follows: (1) shallow soil layer (0—60 cm): the 
isotopic composition and soil water content were highly unstable with time and depth due to 
precipitation inputs and evaporation; (2) middle soil layer (60—120 cm): variations of the isotopic 
composition and soil water content were milder than those in the shallow soil layer; (3) deep soil 
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layer (120-200 cm): the isotopic composition and soil water content were relatively stable over 
time and depth. Source increment was defined as 1% and mass balance tolerance was defined as 
0.1% . Our calculated results are reported as the distribution of feasible solutions and the mean on 
each sampling date. 

Statistical analysis was performed using SPSS 19.0 software (SPSS Inc., Chicago, USA). 
One-way ANOVA was used to assess the temporal patterns of 5!8O values and soil water contents 
in different soil layers by analyzing each sampling plot individually across the sampling periods 
at each soil depth (P<0.05). 


3 Results 


3.1 Precipitation distribution and isotopic variations 


The total precipitation was 86.9 mm from October 2014 to September 2015 (Fig. 3). Less than 5.0 
mm of rain fell during most (91%) of the precipitation events. Precipitation mainly occurred in 
April, July, and September. There were three heavy precipitation events with the amount >10.0 
mm occurred during the sampling period. Specifically, the first one occurred on 17 January 2015, 
with a snowfall of 10.2 mm (8D= —256.51%o, 5'!8O= —32.88%0; n=13); the second one occurred on 
1 April 2015, with a rainfall of 10.8 mm; and the third one occurred on 8 July 2015, with a 
rainfall of 14.2 mm (8D= —27.76%o, 5!8O= —6.77%o; n=13). The 5D and 5!8O values experienced 
large fluctuations during the sampling period. The ôD values of rainfall ranged from —63.39%bo to 
—15.04%o (with the mean value of —31.33%o), and the 5!8O values ranged from —8.0%o to —3.34%o 
(with the mean value of —6.42%bo.) 
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Fig.3 Temporal variations in precipitation and isotopic composition from October 2014 to September 2015 


3.2 Soil water content distribution and isotopic variations 


As shown in Figure 4a, soil water content changed with sampling time and soil depth. In the 
shallow soil layer (0—60 cm), soil water content was highest in July (2.33%) and then dropped 
dramatically to 1.27% and 1.40% in August and September, respectively. Soil water content 
changed little during May, August and September. In the middle soil layer (60-120 cm), soil 
water content displayed relatively mild monthly fluctuations, decreasing from 1.22% in July to 
0.91% in September. Soil water content was relatively stable in the deep soil layer (120—200 cm), 
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varying from 0.92% to 1.16%. Generally speaking, soil water content in the shallow layer varied 
markedly and was significantly different with those in the middle and deep soil layers (P<0.05). 
The mean soil water content decreased with depth and the values were 1.63%, 1.07%, and 1.02% 
for the shallow, middle, and deep soil layers, respectively. 
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Fig. 4 Mean soil water content (a) and 5'8O (b) in the vertical soil profile (0-200 cm) in different sampling 
months in 2015. Error bars represent standard errors (n=4). 


The 5'8O values of soil water showed obvious variations with sampling time and soil depth 
(Fig. 4b). In the shallow soil layer, the 5'8O values of soil water were lowest in August (—3.13%o) 
and highest in May (3.72%o). At depths of 60-120 cm (middle soil layer), the 5!8O values 
fluctuated moderately over time, increasing from —0.84%o to 3.18%o. The 5!8O was relatively 
stable at depths of 120—200 cm (deep soil layer), with the values ranging from —0.95%bo to 2.61%o. 
Generally, the 5!8O values of soil water did not differ significantly (P=0.097) among different 
soil layers. The mean 5!8O decreased with soil depth and the values were 1.43%o, 0.99%o, and 
0.78% for the shallow, middle, and deep soil layers, respectively. 


3.3 Isotopic composition of plant xylem water and water source analysis 


The 5!8O values of xylem water for T. laxa were lowest in September (—1.99%o; Fig. 5). Xylem 
water was enriched in 5!8O in July (-0.39%o). The 5!8O values of xylem water in July differed 
significantly as compared to those in May, August and September (P<0.05). 

As shown in Figure 6, the local meteoric water line (LMWL) was expressed as: 
5D=6.905!80+5.43 (R?=0.84, n=18). The slope and intercept of the LMWL were lower than 
those of the global meteoric water line (GMWL; 5D=865!°O+10; Graig, 1961). However, most 
5!80 values of rainfall samples were distributed around the GMWL, indicating a slight 
evaporative enrichment during precipitation. Moreover, the 6'°O values of condensate water 
samples were distributed on the GMWL, indicating that the condensate water samples had the 
same water vapor source as precipitation and that no evaporative enrichment occurred during 
condensation (8D= -67.55%o, 8!8O= -—9.43%o). However, the isotopic compositions of 
groundwater, soil water and lake water appeared to be subject to severe evaporative enrichment, 
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all of which were below the LMWL with a regression line of 5D=3.915!8O-31.69 (R?=0.87, 
n=185). The mean value of 5!8O was highest in lake water (6.91%o), followed by soil water 
(0.61%o) and groundwater (—5.37%b). 
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Fig. 5 Variations in 5!80 of plant xylem water in different sampling months. Error bars represent standard errors 
(n=4). Different lowercase letters represent significant differences in 5'*O among different sampling months at 


P<0.05 level. 
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Fig. 6 Linear regression relationships between 5D and 5!8O of rainfall, soil water, xylem water, lake water, 
groundwater, condensate water, atmospheric water vapor and snow water relative to the local meteoric water line 
(LMWL), global meteoric water line (GMWL) and evaporation line (EL) 


From Figure 6 we can see that the 6!8O values of xylem water (mean 5!°O of —1.46%o) were 
distributed close to the 5'8O values of soil water (mean 8!8O of 1.02%o), while they were 
significantly different from those of lake water (mean 5!8O of 6.91%o), snow water (mean 5!°O of 
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—32.88%o), atmospheric water vapor (mean 5!8O of —43.28%o), condensate water (mean 5!°O of 
—9.43%o) and groundwater (mean 6!8O of —5.37%o), indicating that T. laxa mainly absorb water 
from the soil. 

The uptake fraction of soil water from the three potential soil water sources (depths of 0—60, 
60-120 and 120-200 cm, corresponding to the shallow, middle, and deep soil layers, respectively) 
varied during the sampling period (Fig. 7). In May, the mean uptake fraction of soil water from 
the middle and deep soil layers was 46.70% and 43.96%, respectively, while it was only 9.35% 
for the shallow soil layer. However, in July and August, the mean uptake fraction of soil water 
from the shallow soil layer increased dramatically while that from the middle and deep soil layers 
decreased, compared with the situation in May. The mean uptake fraction of soil water from the 
shallow soil layer was 44.50% and 40.90% in July and August, respectively. During September, 
the deep soil layer was the main water source for soil (43.33%). 
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Fig. 7 Variations in uptake fraction of soil water from the three potential soil water sources (depths of 0—60, 
60-120 and 120-200 cm, corresponding to the shallow, middle and deep soil layers, respectively) in different 
sampling months. Error bars represent standard errors (n=4). 


4 Discussion 


4.1 Variations in isotopic compositions of precipitation and soil water source 


The differences in the isotopic composition of precipitation may result from the varied water 
vapor origins, atmospheric flow paths of vapor trajectories, temperatures, and precipitation 
amounts (Dansgaard, 1964; Tian et al., 2001; Liu et al., 2014). In our study area, the '8O values 
of precipitation were closely linked with local precipitation amount, especially during heavy 
precipitation events. The 6'8O of precipitation was enriched in June and August when 
precipitation was scarce, while it was depleted in July and September when precipitation was 
abundant (Fig. 3). This result is consistent with the findings of Wu et al. (2016). 

Soil water content in the shallow soil layer (0—60 cm) varied greatly, indicating the influence of 
precipitation recharge, especially the heavy precipitation events, on soil moisture (Allison, 1982). 
The influence of a heavy precipitation event on soil water content was generally greatest in the 
shallow soil layer and was decreased with soil depth. As shown in Figure 4, soil water content in 
the middle and deep soil layers (below 60 cm depth) showed no significant (P>0.05) monthly 
variation, indicating that precipitation had little influence on soil moisture in these layers. The 
5!8O values of soil water also showed marked variations in the shallow soil layer and were more 
positive than those in the middle and deep soil layers due to strong evaporation (Gazis and Feng, 
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2004; Li and Zhou, 2004). This suggests that the shallow soil layer exhibited positive isotopic 
values when soil water content was elevated. However, the middle soil layer displayed less 
variations in soil water content and isotopic composition than the shallow soil layer (Fig. 4), as 
the influence of soil evaporation and precipitation infiltration gradually diminished with soil 
depth. The variations of soil water content and isotopic composition were closely linked with the 
amount of precipitation, antecedent soil water, and soil evaporation. In July, soil water content in 
the shallow soil layer increased dramatically due to a heavy rainfall event (14.2 mm) on 8 July. 
However, the isotopic composition was more positive in the shallow soil layer because of the 
effect of strong evaporation (Ma et al., 2012). 

In desert environments, precipitation is extremely scarce and cannot form runoff, and soil water 
is mainly recharged by precipitation and groundwater. Figure 6 shows that the 6D and 5!8O values 
of some soil water samples were distributed close to the LMWL, indicating that precipitation 
could recharge soil water on occasions. However, the depth of precipitation infiltration was 
limited, and the maximum depth cannot exceed 50 cm (Yang et al., 2014). Therefore, precipitation 
can be a water source only for the shallow soil layer. Groundwater can recharge the deep soil 
layer due to capillary action. Generally, the height of capillary water rise is very limited in the 
desert (<60 cm) because of the presence of coarse soil particles (Zhao et al., 2010). However, the 
5D and 5!80 values of soil water and groundwater clustered along the evaporation line (EL; Fig. 
6). Moreover, the isotopic values of groundwater (5D= —55.63%o, 8!80= —5.37%o) were lower than 
those of soil water (D= —28.96%o, 8'8O=0.61%o). This suggests that soil water may be 
replenished by groundwater in the form of gaseous water or film water, and that the isotopic 
values of groundwater are enriched during the water vaporization process (Chen, 2004). 


4.2 Water use patterns of Tamarix laxa 


The isotopic composition of xylem water exhibited seasonal variations, suggesting that the water 
absorbed by T. laxa species from May to September originates from different soil layers (Fig. 7). 
In July, the xylem water had less negative isotopic values as compared with the soil water in the 
middle and deep soil layers, while these values were similar to those in the shallow soil layer. 
This implies that T. laxa species mainly absorbs water from the shallow soil layer after heavy 
precipitation events. During the study period, a heavy precipitation event (14.2 mm on 8 July) 
occurred before sampling (Fig. 3), and soil water content largely increased in the shallow soil 
layer (mean value increasing from 1.53% to 2.33%) while exhibited no apparent variations in the 
middle and deep soil layers (Fig. 4). In addition, T. laxa species has a developed root system in 
the shallow soil layer (Fig. 2). Therefore, we conclude that T. laxa mainly absorbs water from the 
shallow soil layer. In contrast, xylem water of 7: laxa showed more negative isotopic values 
during other months (May, August and September). At the beginning of the growing season, 
xylem water displayed more negative isotopic values, similar to soil water in the middle and deep 
soil layers. A heavy precipitation event (10.8 mm) occurred on 1 April, during which soil water 
content in the shallow soil layer was replenished; however, this event was not as striking as that 
(14.2 mm) in July and soil water content was far from the field water capacity. Additionally, due 
to the low temperature and water consumption at the beginning of the growing season, T. laxa did 
not utilize water from the shallow soil layer and instead depended on water from the middle and 
deep soil layers. There was little precipitation from mid-July to the end of August. Moreover, 
evaporation was severe during this period. Thus, soil water content of the overall soil profile 
decreased markedly. However, the shallow soil layer could remain wet for a short period after an 
exceptional precipitation event, because of a dry surface sand layer preventing evaporation loss 
(Liu et al., 2006). Consequently, T. laxa could continue to utilize water from the shallow soil layer 
in August. Although there were many precipitation events in September, the recorded value for 
each event was less than 5.0 mm (except for one on 3 September, with a recorded value of 5.5 
mm), which is generally not sufficient to infiltrate into the soil. Consequently, soil water content 
in the shallow and middle soil layers did not demonstrate any recharge at this time, but instead 
decreased gradually due to strong evaporation. During this stage, T. laxa mainly absorbed soil 
water from the deep soil layer. 
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5 Conclusions 


Stable isotopic techniques were used to explore the water sources and water use patterns of T. 
laxa in the mega-dunes of the Badain Jaran Desert. Due to strong evaporation, the isotopic 
composition of soil water was more enrichment in the shallow soil layer (0—60 cm) than in the 
deep soil layer (120—200 cm). Precipitation amounts, especially from heavy precipitation events, 
strongly affected the vertical distribution of soil water content and isotopic composition. We 
found that the values of isotopic composition in snow water, atmospheric water vapor and 
groundwater were more negative, while those in lake water were more positive. The isotopic 
compositions of all these water sources differed from the isotopic composition of xylem water in 
T. laxa. The major water source of T. laxa is soil water, although soil water content varies 
throughout the year along the depth. The water use patterns of T. laxa reflect its ability to absorb 
water from different soil layers according to soil water content. Precipitation could infiltrate into 
the shallow soil layer but cannot contribute to groundwater recharge. 
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